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We develop a microscopic theory of sound damping due to Landau mechanism in dilute gases 
with a Bose condensate. It is based on the coupled equations that take into account the mutual in- 
fluence of condensate density, phase (or superfluid velocity) , and distribution function of elementary 
excitations on evolution of the system. These equations have been derived in earlier works within 
a microscopic approach which employs the Peletminskii-Yatsenko reduced description method for 
quantum many-particle systems and Bogoliubov model for a weakly nonideal Bose gas with a sep- 
arated condensate. The dispersion equations for sound oscillations are obtained by linearization 
of the mentioned evolution equations in the collisionless approximation. They are analyzed both 
analytically and numerically. The expressions for speed of sound and decrement factor of sound are 
obtained in the limiting cases of high and low temperatures. It is shown that at low temperatures, 
the decrement of collisionless damping of sound has a quadratic dependence on temperature. Such 
dependence essentially differs from that obtained earlier by other authors in phenomenological ap- 
proaches. It is also demonstrated that at high temperatures, the sound damping coefficient has a 
linear dependence on temperature that coincides with earlier results of other authors. We discover 
a finite temperature at which the speed of sound is exactly equal to the speed of sound at zero tem- 
perature. This effect is due to non-analytic dependence of dispersion characteristics of the system 
on temperature at this point. We also indicate the possibility of oscillations for the decrement or 
increment coefficients of sound in the vicinity of the discovered temperature. 
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I. INTRODUCTION 

The study of mechanisms of sound damping in a Bose- 
Einstein condensate (BEC) has a long history. Calcu- 
lation of the sound damping coefficient in systems with 
BEC is a rather complicated theoretical problem. First 
expressions for damping rate in such systems have been 
apparently obtained in [1 Q for the spatially homoge- 
neous case. 

The direct experimental observation of BEC has 
simulated a great number of works devoted to various as- 
pects of this phenomenon (see in this regard, for example, 
0,0] and references therein). A number of papers, both 
theoretical and experimental, deal with the problem of 
propagation and damping of excitations in Bose gases 
with the presence of a condensate [8l4l6|| . 

It is currently assumed that the Landau damping is the 
most probable mechanism of sound relaxation in the so- 
called trapped Bose condensates. This mechanism con- 
sists in collisionless absorption of the energy of oscilla- 
tions by quanta of elementary excitations |12l . [l3j . In 
this regard it is worth recalling that the existence of spe- 
cific collective excitations in a gas with BEC is known 
since the pioneering work of Bogoliubov [l7j . In this pa- 



per, a special perturbation theory has been proposed for 
a weakly non-ideal and spatially homogeneous Bose gas 
with a condensate in which the repulsive interaction acts 
between atoms. This theory predicts the elementary ex- 
citation spectrum for such system at zero temperature. 
At small wave vectors, it coincides with the spectrum of 
sound oscillations in a condensed Bose gas. 

The first work, which proposed a method for calcu- 
lating the sound damping rate in trapped BEC due to 
Landau mechanism, is apparently a paper [l2j (see also 
[l3j). The approach developed by the authors uses the 
perturbation theory and is based on the calculation of the 
difference in probabilities between emission and absorp- 
tion of quanta of oscillations by elementary excitations in 
the system. It should be noted that the trapped conden- 
sates represent a spatially inhomogeneous system. This 
fact essentially com plic ates the analytical calculations. 
The method of Ref. [12| has been shown to be suitable 
for numerical calculations of sound damping in a trapped 
condensate (see the same paper). For a spatially homoge- 
neous BEC, this method has provided analytical formulae 
for damping according to Landau mechanism (l2j . In this 
case the authors have reproduced the results obtained in 
Refs. [10. 
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It should be stressed that the formulae of Ref. [121 ] 
(and, hence, the results of Refs.JT], [1]) can be obtained 
by other method involving the kinetic equation for dis- 
tribution function of elementary excitations (in this con- 
nection see, e.g., Ref. [6j). In other words, the semi- 
phenomenological approach to calculation of the damp- 
ing coefficient shown in Ref. [f| is equivalent to the 
method of Ref. |l2|. This is inter alia indicated in 
Ref. [Hj]. Thus, the correctness of the problem of sound 
damping in BEC which is under solution becomes de- 
pendent on the type of kinetic equation for elementary 
excitations. Therefore, we must be confident of correct- 
ness of kinetic equation which is used to solve the men- 
tioned problem. However, it is clear that in the most 
general case the only use of kinetic equation for excita- 
tions is not sufficient. The system should be described 
by the coupled evolution equations which take into ac- 
count the mutual influence of condensate density, phase 
(or superfluid velocity) and distribution function of ele- 
mentary excitations. Thus, there is a necessity for con- 
sistent (controlled) derivation of the coupled equations 
describing the system with BEC at a particular stage 
of its evolution. The consistent derivation of the sys- 
tem of coupled equations can be achieved using a mi- 
croscopic approach, proceeding from the first principles. 
This problem has been solved in [l(J H 01, El within 
a microscopic approach based on the reduced descrip- 
tion of quantum many-particle systems and Bogoliubov 
model for a weakly non-ideal Bose gas with condensate 
17]. Note that the foundation of the reduced description 
method for classical (not quantum) systems has been for- 
mulated by N. Bogolyubov 22j. This method has been 
extended to quantum many-particle systems by S. Pelet- 
minskii and coauthors (see Ref. 1231 V The synthesis of 
approaches elaborated in Refs. \l7\ and [23[ has allowed 
one to obtain the kinetic equation for distribution func- 
tion of elementary excitations coupled with the evolution 
equations for condensate density and superfluid momen- 
tum [l8|, Oil ■ The validity of such system of equations is 
confirmed by its controlled derivation within the frame- 
work of special perturbation theory in weak interparticle 
interaction. In addition, the following fact says in fa- 
vor of the mentioned coupled equations: they have been 
employed to derive hydrodynamic equations of a super- 
fluid in which the smallness of the difference between su- 
perfluid and normal velocities is not taken into account 
pM |2l| . As this difference tends to zero, the obtained 
equations are reduced to the well-known hydrodynamic 
equations of Khalatnikov (see e.g. [24[). 

Notwithstanding the above considerations, the equa- 
tions of Ref. [19j have not yet been used to study the 
propagation and damping of sound in a dilute gas with 
BEC. However, the same considerations allow us to hope 
that the correct solution of evolution equations found in 
Ref. [l9| should lead us to the correct expression for 
the the sound damping coefficient in a gas with BEC. 
The present paper is devoted to the study of collisionless 
mechanism of sound damping in dilute gases with BEC 



on the basis of general dynamic equations of such sys- 
tems obtained in [19j from the first principles. As will be 
seen later, the results obtained in the present work signif- 
icantly differ in some cases from those of Refs. [l|, [|[ and, 
consequently, of Ref. [l2[ • For example, in the present 
study it is shown that in the collisionless approximation 
the sound damping factor at low temperature is propor- 
tional to the square of the temperature 7 ~ T 2 , whereas 
the results of |l|, [H| give 7 ~ T 4 dependence. 



II. THE KINETICS OF SPATIALLY 
INHOMOGENEOUS BOSE GAS IN THE 
PRESENCE OF CONDENSATE 

In the construction of the kinetics of spatially inho- 
mogeneous Bose gas the authors [l9[ started from the 
Liouville equation for the statistical operator p{t) 



. dp jt) 

1 at 



H,p(t) 



(i) 



where H = H + V is the Hamiltonian of the system, 
consisting of the ideal gas Hamiltonian Hq 



H = — [d 3 r d ^ +{r) d ^ [r) 
2m J dr k dr k 



dr k dr k 
and binary interaction Hamiltonian V 



(2) 



V = - J d 3 r J d 3 RtP + (r + R) (r) V (\R\) g 
x-0 (r) i)j (r + R) . 

Equation (fT]) is written in units in which Planck's con- 
stant h is equal to unity. In the formulae ([2]), (|3]) m is 
boson mass, V^d-RI) is the binary interaction potential, 
which depends only on the distance between particle, and 
(r) , ip (r) are field operators. As the parameters de- 
scribing the weakly non-ideal Bose gas with condensate 
in the kinetic stage of evolution of the system in [19], the 
quasiparticle distribution function f p (r,t), the order pa- 
rameter %l>(r,t) — Spp(t)ip (r) and superfluid velocity 

Vk ( r ,t) = m^ 1 -£^-Imln Spp(t)iil! (r) were selected. 

Using the method of reduced description [23[ combined 
with the special perturbation theory [I?} allowed to ob- 
tain in [191 ] the following system of equations for the pa- 
rameters f p (r, t), ip (r, t) v k (r,t): 



dfp 
dt 



+ Q P (f,i>)=L p 



dip 



dtp I dv k 
dt ' Wfc ^ + 2%^ =M/>) 



(4) 
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where the value L p (/, ip) is the quasiparticle collision 
integral and functionals (f,tp), Q p (f,ip), h(f,ip) are 
given by the expressions: 



W (/, ip) 



d 3 p 



Or 



2ip J (2tt) j ^ 



L p (f,iJj)+o(X 5 ), 



df d df d 



dr k dp k 

dbJr, 



1 



2mip \ dip 



dpk dr k 
va I a 
2 \ x 



A+h (V^)+o(A 4 ), 
p I 

(5) 



L p (f°,if>)=Q, i^(/°,V)=0, 



u p — pv 



T 



(8) 



with the chemical potential of quasiparticles is equal to 
zero; here T is the temperature in energy units. The 
vanishing of the chemical potential reflects the fact that 
the number of the quasiparticles is not conserved during 
collisions 0,[2l|. 

In the collisionless approximation equations (HJ) - © 
can be written as 



where o(X n ) is the value in order of magnitude A™ and 
A characterizes the smallness of the interaction between 
the particles, ~ A 2 . Let us recall that the basis 

of Bogoliubov equilibrium state theory of Bose gas in the 
presence of interaction is the assumption that ip ~ A , 
see [13] • Furthermore, it is believed that the order of 
magnitude of ip (r, t) does not change after differentiation 
ip(r,t) with respect to r, (dip (r, t) /drk) ~ A -1 . In 
the formulae (HI ([5]) the following notations were also 
introduced: 



dfp 
dt 



df p d dfp d 

or k dp k y dp k dr k 



dn 
~di 



9 , x 
dr k 

^r + dr- k \Y + hif ' n < 



(9) 



0, 



where we have introduced the new variable n(r,t) that 
is the condensate density [2l| 



Bp (ip,v) = u p (i/j)+pv, 
(A p ) = J d 3 pf p Ap, u = f p=0 , 

vp = J d 3 RV (\R\) exp (-ipR) , 



(6) 



p ^p 



p#0 



v^ip 2 



where A p is an arbitrary function of momentum p. 
Appearing in ([5]) and (j6|) values are given by formulae: 



up ty) = [a 2 p (iP) - I3 2 p (J,)] 
a p (ip) = e p + P p (tp) , 



1/2 



Zm m 



(7) 



We note that quantity u p represents Bogoliubov spec- 
trum of elementary excitations. 

The explicit form of the quasiparticle collision inte- 
gral is not provided in the present paper because it is 
not necessary. As mentioned earlier, we will explore the 
mechanism of collisionless sound damping in a gas with 
BEC (i.e. Landau mechanism). The explicit form of the 
quasiparticle collision integral can be found in [lj| . 
Here we note only the fact that the quasiparticle collision 
integral L p (/, ip) (and hence (f, ip), see (J5])) vanishes 
by substituting stationary Bose distribution function /° 



ip 2 (r,t) ee n(r,t), 



(10) 



and the values e p (n, v) and h (/, n) are given by expres- 
sions 



£ p (n, v) ee uj p (n) + pv, 

h (f,n) = I \ + 21 /<) + hQ (n) ; 

m \ u n Zm \ u„ 



J p I 

e p - Up ap - Up 
Vp V Vq 



m ' 
(11) 



where 



Up (n) ee [e p (s p + 2/3 p (n))] 1/2 , 

u n (12) 

a p (n) ee e p + p3p (n) , pJ p (n) = , 

m 

and the value e p is still defined by ([7]). Equations © to- 
gether with @ and (fTT|) (fT2|) and are the initial system 
of equations to study the propagation and damping of 
sound in dilute gas with condensate neglecting the colli- 
sions of the quasiparticles. 



III. SOUND DISPERSION EQUATIONS IN 
DILUTED GAS WITH BEC 



In order to investigate the propagation of the sound in 
jas with BEC, we linearize the system of equations ©, 
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(fTTjl and (|T2]) with respect to a spatially homogeneous 
equilibrium state according to following formulae 

n(r,t) = n + h(r,t), 

v(r,t) = v(r,t), (13) 
f P (r,t) = f° p +f P (r,t), 



n > \n{r,t)\, f° > / p (r,t) 



where no = no (T) is the equilibrium value of the con- 
densate density in the system at the temperature T. The 
equilibrium value of the velocity vo is considered to be 
equal zero in the second formula in (1131) . Thus the ve- 
locity |S(r,t)| is supposed to be in order of magnitude 
of h(r,t) and f p (r,t). Then the equilibrium distribu- 
tion function of quasiparticles pursuant to (jU) is given by 
expression 



/« = (cxp^-l 



(14) 



In this formula as in all subsequent expressions we 
omit the index "0" in the designation of the equilibrium 
value of ujp = lu p (no) to avoid encumbering of com- 
putations. Note that the quantity h (/ , no) defined by 
(ITT1) represents the chemical potential of atomic Bose gas 

= h (f°,n ) (see in this regard 

Deviations of the corresponding quantities from their 
equilibrium values were denoted by fi(r,t), v(r,t) and 
f P (r,t). 

The equations of motion for these variables can be rep- 
resented as follows: 



Proceeding further to the Fourier transform of quan- 
tities h(r,t), v(r,t) , f p (r,t). in equations (|T5|) as pro- 
vided by formula 



diu / dV 9r - iwt C(<Z^), (17) 



where ((r,t) should be understood as either of the con- 
sidered quantities, we obtain 



um(q,Lu) = n a qv (q,u) , 



Vq 



uv k (q,uj) = —q k n(q,uj) + q k A (q, w) , 



W - q k 



du} p \ 
dp k ) 



(18) 



-Qk 



dp k 



pv (q,uj) + ^^-n(q,u}) 



where we denote (see (fl~5|) . (fT6| ) 



1 



A(q,u) = J d 3 pf p (q,u})K p 

e p (vo + 2v p ) + v p vono 



2m (2tt) 



3 / d 3 pf p (q,u) 



(19) 



In the formulae (fT5)l (fT9"f in the Fourier transforms of 
the quantities n(r,t), v(r,t) , f p (r,t) we omit the 'tilde' 
sign. 

Further expressing the values n(q,uj) and v (q,uj) from 
the first two equations of (fT8|) in terms of A (q, co), 



d . d _ , . 

—n(r,t) + n — v fc (r,i) = 0, 



d . vo d _ 

—v k (r,t) H — - 

at m arfe 



ft(r,t) + /" d 3 pK p ^-f p (r,t) = 0, 



where we have introduced the following notation 
1 e p (vo + 2v p ) + v p v no 



K„ 



'" 2m (2tt) 3 



(15) 



(16) 



Recall that in the present paper Planck constant h is 
considered be equal to unity. Note also that deriving 
equations (|15[) we have discarded terms that are higher 
than the first order of magnitude A. It was taken into 
account the circumstance of (JTUJ) that n ~ A~ 2 or i.e. 
v p no ~ A . 



( \ n °1 A I \ 

n(q,u) = -5 Y^A(q,uj) , 

or — Mq<7^ 

I \ W * A I \ 

v k (q,uj) = — T ^A(q,uj) , 



(20) 



the third equation of (|18[) can be written in the form 



: p no 

uj p (cj 2 - ulq 2 ) 



dp k dw p 



pqujujp + v p £ p noq" . . . 



(21) 



where u is so-called speed of zero sound in Bose gas 
2 _ v n 



(22) 



As is readily seen, equation ([211) subject to (fl9| is 
an integral equation for distribution function Fourier 
transform/^ (q, ui). The solution of this equation can be 
represented as in [25| : 
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fp («,<*>) = S p (q)8(uj - qu p ) 



qu 



df° v 



pquJUp + v p e p n q 



T^A(q,Uj), 



p dujp uip (uj 2 — Mq<7 2 ) (uj — qu p + iO) 

where we have introduced the following notation: 

_ dcjp 
Up = ~dp~' 



(23) 



(24) 



df° 

Ia V (q, w) = B ap (q) 5(u>- qu p ) + e 1 (q, ui) Q u Pq^- 

pqujujp + VpEpHoq 2 
uj p (u z — uffi) (uj — qu p + iQ) 

(29) 

We remark that in case of charged particles gas the 
value e{q,uj) (see ([2"5])) 



and Bp (q) is an arbitrary function, which is required to 
impose the following restriction: the distribution func- 
tion f p (r, t), calculated in accordance with flT3"l) . (fTT)) and 
(|23p should be small compared with the equilibrium dis- 
tribution function /°. 

It implies also from (IT7|) that B p (q) must satisfy the 
following relation: 



e(q,uj) =e 1 (q,u})+ie^(q,u)) 



(30) 



represents complex dielectric permittivity of the system 
(see e.g. [23j]). It is known that the presence of imaginary 
term in dielectric permittivity indicates the energy dissi- 
pation of electromagnetic waves with dispersion relation 
that should be obtained from the equation 



B;(q) = Bp(-q) 



(25) 



e(q,uj a (q) - i lq ) = 0. 



(31) 



All the valid set of such functions we denote by 
B a p (q), where a is a continuous or discrete symbolic pa- 
rameter such that the functions B p (q) = B ap (q) may 
depend on a. The reason for introducing this index may 
consist, for example, the following: the set of functions 
B av (q) should be enough to build an arbitrary value for 
the distribution function f v (r, t) at the initial time t = 
(see in this context also |26|). 

Formula (f2"3")l allows to find the value A (q, uj) in terms 
of B av (q) functions: 



Aa (q,u) = Ba (q,uj)e 1 (q,u) . 



(26) 



where 



Moreover, the wave decrement j q is determined by the 
imaginary part £2 (<Z,w) of value s (q,uj). For this reason, 
weakly damped oscillations in the system 



N (<?)| > lq 



can exist if only 



|ei (q,ui)\ > \e 2 (q,u)\ 



(32) 



(33) 



besides, as the consequence of (f30|) - ([33]) (see in this 
regard [25, 26]), the frequency ujq (q) can be found from 
the equation 



Ba (q,u) 



3 VQflQVp + e p (v Q + 2v p ) 
3 , a p 



2m (2ir) 

xB ap (q)S(uj - qu p ) 



(27) 



and 



e(q,u) = l-\ ^—-3 / d 3 p—^- (qpujuip + v p E p n q 2 

2m (2tt) J clujp 

qu p [v noVp + e p (vq + 2v p )} 



ujp (uj 2 - Mg<7 2 ) (uj - qu p + iO) ' 



(28) 



The expression (|2"3"j) for given values f p (q,uj) subject 
to (|26p - (|2"81 can now be represented in the form: 



ei(g,w (g))=0 (34) 
and the damping rate 7 9 is given by expression 



dei (g,ui) 
dui 



e 2 (q,uj (q)) . 



(35) 



U— U)Q 



The structure of solution (|29|) is such that weakly 
damped waves may also be excited in investigated sys- 
tem, in accordance with (f2"U|) and (fT?) . and the disper- 
sion law is determined by (|3"Tj) - (f3"3")l . As we will show in 
the following section, such waves would represent sound 
waves in weakly nonideal Bose gas with condensate. The 
value 7 9 obtained according to ([29)) - (|32j) will determine 
the attenuation of sound in investigated system. 
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IV. THE DAMPING RATE OF SOUND IN 
DILUTE BOSE GAS WITH CONDENSATE 



After the same angle integration the second of expres- 
sions (f3"Tf can be written as follows 



In order to solve this problem it is necessary to deter- 
mine the real £\ (q,uj) and imaginary e 2 (q,uj) parts of 
the value e (q, uj). For this purpose in the expression (j28|) 
we use the formula 



= P in5 (x) , 

x + iO x w 



(36) 



where the symbol P means that the further integration 
is taken in means of Cauchy principal value. After some 
hackneyed transformations we obtain the following ex- 
pressions for £1 (q, w)and £2 (q, uj): 



£1 (q,u) = 1 



1 



2m (2tt 



,p / 



duj v 



qu p (qpujuip + f p £ p n g 2 ) \vqtiqv p + e p (V + 2v p )] 
uj 2 (uj 2 - ulq 2 ) (ui - qu p ) 



p 

71 id 



jj f Of 

£ 2 (q,uj) = - " / d 3 p—^ (qpujuj p + v p e p n q 2 ) 

2m (2n) J ouip 

[u n v p + e p (isq + 2v p )} 



uij (uj 2 - u 2 q 2 ) 



-5 (uj — qup) 



(37) 



Further calculations is not possible without specifying 
the explicit form of v p that is the Fourier transform of the 
interaction potential y(|iZ|) (see ([3]), ©). To simplify 
the calculations, in many cases it is accepted to replace 
the interaction potential V ( \ R]) by the following effective 
potential (see, in this regard 



6. 7]) 



£2 (q,u) 



v uj 



00 dpp 2 dfp 



87rm (uj 2 — u^q 2 ) J Q ui 2 u p dui p 

x [(p/u p )uj 2 ujp + e p v n q 2 ] (n v + 3e p ) 9 (qu p - \uj\) , 

(41) 

where 9 (x) is Heaviside step function. 

We now define, following the work [ijj], the sound 
damping rate in the case of low (T <C vqUq) and high 
(T ^> vqtiq but T < T c ) temperatures. To study the 
case of low temperatures it is convenient to use variable 
z = uip/Tin the integration in expressions (|4"0")) and (|4T|) . 
The variables p, e p , u p in terms of the variable z can be 
expressed using the u> p explicit expression (|12[) and also 
([7]) and (j3"9")l . Carrying out all rather cumbersome but 
necessary calculations, the result can be summarized as 
follows: 



£1 (q,uj) = 1 

£2 (q,u) 



F 



a(uj 2 ~u 2 q 2 ) JQ 



dz — g(q,uj,z) 

oz 



2quov (z) 
ttF 



In 



uj + qu v (z) 



ui — quov (z) 



, df° (g) , s 
2a(uo 2 -uW)h dz ^T 9{w) 



qu v (z) 



e(qu v(z)-\uj\),f (z) 



(42) 



where the functions g (q,uj, z) and v (z) are given by: 



V(\R\)=v 5(R), u Q 



Ana s 



(38) 



where a sc is so-called scattering length (for details see 
[6]). Thus we have v p = vq. Taking this into account we 
have ujp = uj p and the value u p (see (1.12) (|2"4"|) ) can be 
expressed in the form 



Up — Wp, 

v 



dujp 
dp 



g p + vonp p 
uj v m 



(39) 



Formulae (1381) and (|39j) enable us to perform integra- 
tion in ([37]) over the angle between vectors q and p , and 
then £1 (q, uj) can be represented in the form 



£1 (q,uj) = 1 
V 



dpp 2 df° 



mu 



v 

xb 



An 2 (uj 2 - u^q 2 ) J uj 2 8uj p 
(n is + 3e p ) 
u + qu p 



UJ 2 UJ P + e p u 2 q 2 



(40) 



2qUp 



In 



qu p 



<(q,u),z) 



V2VVTTW -1 (3Vl + a 2 z 2 - 2) 



z (1 + a 2 z 2 ) 



q 2 u, 



2 (l + a 2 z 2 - Vl + a 2 z 2 ) 



222 
uj a z 



v(z) 



V2Vl + a 2 z 2 VVl + o?z 2 ~ 



az 



(43) 



and the following notations were introduced (recall uq 
still given by 



(44) 



Deriving the last formula of (|44[) it is necessary to use 
expressions (|22p and (I38J) . The quantity Ho a L i s called 
the gas parameter, seejg, 0|- That is, in diluted gases 
the following relation should be satisfied: 



T 

a = , 

v n 

I 2 
F = —-v m 2 u = —=yjn al c . 
47T- 2 Jit 
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n a sc <C 1. 



(45) 



It is easy to verify that the values E\ (q, oj) and £2 (g, u>) 
can be represented as functions of one dimensionless vari- 
able s = lu (q) /uoq 



ei (s) = 1 



F 



a (s 2 



1) 



1 - 



£2 (s) = - 



2v (z, a) 
irFs 



hi 



' df° (z) 
dz—^—g{s,z) 

■f ii (z, a) 



2a (s 2 



1 



s — v (z, a) 

df°(z)g(s,z) 



(46) 



dz- 



dz v(z,a) 



x6 (v (z, a) — \s\) 



where 



g(s,z) 



2VV1 + a 2 z 2 - 1 (3\/l + a 2 z 2 - 2) 
z(l + a 2 z 2 ) 



(47) 



1 + a 2 z 2 - \Jl + a 2 z 2 + s 2 a 2 " 2 



and for these functions the symmetry conditions are valid 



£i(s) = si(-s) , e 2 (s) = e 2 (-a) . 



(48) 



Expressions (|34j) and (|35|) that determine the disper- 
sion and decrement (or increment) of small oscillations 
in investigated system, can be written in the following 
form by taking into account (|46|) : 



I \ n _ w (<?) 



7? = u aq { £2(s) 



ds 



u q 



As is easily seen, these oscillations have linear spec- 
trum. In virtue of (j4"rJ)) and (|49p the structure of the 
dispersion equation is such that the unknown value s, 
which determines the oscillation frequency as a function 
of the wave vector, does not depend on wave vector itself. 

In accordance with and (|4*^)l the dependence of 
the oscillation frequency ujq on the wave vector q in this 
system should be determined by the solution of equation 



- 1 = — / dz 



1 



so 



2v(z[ 



In 



df°(z) 

dz 
s + v (z) 



g(so,z) 



s - v (z) 



(50) 



We note that dispersion equation ([50|) is similar to the 
dispersion equation of zero sound in a normal Fermi liq- 
uid, (compare e.g. with the corresponding formula in 
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0.99985 - 
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FIG. 1: Dependence of dimensionless speed of sound s (a) = 
loo {q,a)/uoq on temperature (dimensionless quantity a, see 
(144)) ') obtained by a numerical solution of (f50|) for F = 10~ 3 
(solid line). Dots display intersection with line s — 1 (dashed 
one) in points a = and a w 0.847. 



25]). Further, taking this equation into account to calcu- 
late the derivative 9e g^ appearing in (|4*5|) . the damping 
rate 7 g can be presented as follows: 



7g = 
u q 

where 



Fso 



2 ab(s ) J 



dz— -— 9(v(z) 

oz v [z) 



kol), 

(51) 



b(so)^\(s 2 -l) 



dei (s) 
ds 



'0 



1 



(49) + £±£ J JL. 



F 

s a 
Fs f d 



— / dz- 



df°(z) (g(s ,z) dg(s ,z) 



dz 



2a 1 dso 



o 

°° d _ df°(z)g(s ,z 
dz v (z) 



In 



so 



ds 
v(z) 



so - v (z) 



(52) 



We emphasize that mentioned above expressions (|42|) 
- ([52^1 are exact, despite the fact that we have modified 
them to the form suitable to study the low temperature 
regime. 

As is readily seen, equation ([5Tfl) in the general case can 
be solved only by numerical methods. Figure 1 shows the 
dependence s (a) (see (|46l) and obtained as a result 
of numerical solution of equation (|5T)|) for F = 10 -3 . It 
is evident that the function s (a) behaves nonmonotoni- 
cally as a changes. At low a region (the case of low tem- 
perature) increase of function is observed. At the point 
a 0.43 it reaches maximum, and then s (a) decreases 
monotonically as aincreases, and in the point a 0.847 
s (a) equals unity again. As mentioned above, the in- 
crease of value a is restricted at least the critical temper- 
ature, see (|31]). At zero temperature (a = 0) we deal with 
a classical zero sound in Bose system, because s = 1 , and 
hence wq (q) = uoq due to (|46p . This result is naturally 
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expected. However, we note one more time that at the 
point a s» 0.847 also holds s = 1 and the frequency of 
sound in dilute gas with Bose condensate equals again to 
the frequency of zero sound in such a system. The fact 
that we discovered is still not mentioned in the literature. 
It may be a consequence of use of the evolution equations 
of the system in present work, that were obtained in the 
microscopic approach, as well as the numerical solution 
of the dispersion equations. Whatever the case is, the 
very existence of the second point of a sound dispersion 
curve with s = 1, that is a w 0.847 (or T 0.847 v n Q ) 
for F = 10 -3 , requires an individual physical interpreta- 
tion. Some more comments regarding this point will be 
given below. 

In the regions of low (T <C von-o) and high (T 3> vono 
but T < T c ) temperatures a solution of (|5U|) can be ex- 
pressed in analytical form. Thus one can obtain analyti- 
cal expressions for the quantity "f q in two limiting cases. 
To do this, as we shall see, one need to involve a numer- 
ical analysis as an auxiliary technique. 

Consider first the case of low temperatures. By virtue 
of inequality (low temperature regime, as mentioned 
above) 



T 



< 1, 



(53) 



and the rapid decrease of function / (z) as z — > oo (see 
(|42]> ). functions g(q,uj,z) and u(z) can be expanded in 
powers series of az in (|4"6")l - (14*51) : 



When F < 1 and a < 1 (see (|H|). ([15])) the solution 
of ([5B]) can exist only in the 5s <C a 2 -C 1 domain. This 
inequality makes it possible to neglect the value 5s in- 
side logarithm term in the integrand in (I56[) as the first 
approximation of perturbation theory. As the result we 
obtain 



(57) 



5s = Fa 2 d, 5s = s — 1, 
where the notation d is 



Formula (|57[) is also confirmed by numerical calcula- 
tions. Fig. 2 shows the dependence d(5s) , 5s = s — 1, 
computed according to the formula (157)1) 



d(5s) — — f dz 

4 J ( e *_i 



/o 

<! i In 2 - 1 - I In 
2 2 



c 3 2 2 

ds — —a z 



(59) 



As can be seen from the plot on Fig. 2, the equality 
d(5s) « d w 4.4 for a = 0.1 holds up to the second 
decimal. As a consequence, in the expression for d that 
is given above, the dependence on 5s can be neglected. 
Then one can derive: 



g (s, z) a 3 2 2 



w(z) 



l + -a 2 z 2 . 
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(54) 



As will become apparent from the subsequent formulae 
such expansion corresponds to the development of pertur- 
bation theory with respect to small parameter a. Taking 
into account expansions (|54[) the dispersion equation (|50|) 
can be written as: 



si - 1 = Fa 2 



dz z 



,df°(z) 



dz 



1 - 



•so 



2v (z) 



In 



,s + 1 + |a 2 z 2 



so 



|a 2 z 2 



(55) 



■so 



_ w (g) 



"o<7 



In the region of small a, i.e. low temperatures, as we 
have seen from the Figure 1, s « 1. That is to say that 
value 5s = s — 1 should be much less than unity, |<5s| <C 1. 
Therefore, the equation (|55j) admits further simplification 



5s, 



-Fa' 



dzz 



2 df° (z) 



l-iln2 + -ln 
2 2 



dz 

3 
— < 

8 



jr "22 
OS a Z 



(56) 



lna- 1.27, 



(60) 



where the second term was calculated numerically and 
deriving the first term the following integral value was 
used: 



5 4.407 - 
> 

P 4.406 - 

3 

O 4.405 - 



o {e z - 1) 



7T 

2 = y 



0.0000 0.00002 0.00004 0.00006 0.00008 0.0001 



6s parameter 



FIG. 2: Dependence d (5s) plotted numerically according to 
for a = 0.1. 
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Thus the solution ([57]) shows the following dependence of 
the sound frequency in dilute gas with BEC on the low 
temperature range: 



w (?) ~ ±«o? (l + Fa 2 d) , a 



v n 



< 1 



(61) 



Analytical dependence s (a) in case of a 1 in ac- 
cordance with ([57]) and ([62]) is shown on Fig. 3 for three 
different values of F. 

Damping rate of sound in BEC calculated according 
to formulae ([51)). (|52l and using (J54J) , ([57]) - ([61]) in the 
first order approximation with respect to a is given by a 
fairly simple expression 



7 9 



w (?) 



7g 
tt ? 



7T 3 „ 2 7T 5 / 2 



—Fa 2 = 



\Aoa 3 



T 



(62) 



If we assume that at low temperatures (see ([55]) )) all 
Bose gas is in the condensate [l2[, no (T) « n (0) = 
ri , the quadratic dependence of the decrement 7 g on 
temperature follows from ([6"T]) . This behavior of sound 
decrement at low temperatures significantly differs from 
that given in pj ll2T |. In these papers the value 7 9 is 
known to depend in the low-temperature regime on the 
temperature as T 4 . This is obviously a consequence of 
our use of the coupled evolution equations derived in the 
microscopic approach from the first principles. 

Fig. 4 shows the dimensionless damping rate 7 g /w (?) 
at low temperatures (low values of a) plotted for F — 
10~ 3 . The solid line in the figure reflects the behavior of 
the dimensionless decrement that follows from analytical 
expression ([62]) . The dotted line shows the dependence 
j q /ujQ (?) on the temperature obtained as a result of nu- 
merical calculations based on formulae ([50]) - (|52l) . It 
is evident that analytical expression ([52"]) coincides with 
good accuracy the damping rate temperature dependence 
in a dilute gas with BEC in a rather wide range of low 
temperatures. 




0.02 0.04 0.06 0.08 0.10 



temperature 



FIG. 3: Dependence 8s (a), 5s = uiq (q, a)/uoq — 1, plotted 
according to analytical expressions (|57|l , (I60|l for F\ = 3 • 10 -4 
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FIG. 4: Dependence of dimensionless damping rate "/q/uoo 
on temperature (dimensionless quantity a, see (|53[0 . Solid 
line represents the analytical approach given by (|62[) and dot- 
dashed line represents numerical calculation in accordance 
with ([50]) - ([52]). Both lines were plotted for F = 1CT 3 . 
Dashed vertical line indicated the region where analytical for- 
mula (|53[) works with good accuracy. 



Here we should make the following remark. The ex- 
perimental data (see [1, 0] and [27T - l3lj ]) show that un- 
der present experimental conditions, the criterion a — 
(T /vqUq) <C 1 is generally not realized. For example, the 
approximate estimates of the parameter a obtained in 
accordance with typical experimental conditions of men- 
tioned studies, show the following circumstance. This pa- 
rameter takes on the value a ~ 2for 85 Rb [30], a sa 23 for 
23 Na [4], a w 58 for : H 27]. The values of the scattering 
lengths in order to perform these estimates were taken 
from 0. The exceptions are the data of (28[ and [3lj ]. 
For 7 Li we have a « 0.5 [28[ and for 133 Cs accordingly 
a « 0.8 [3l[. Therefore, to observe the effect of sound 
attenuation in dilute gases with BEC at low tempera- 
tures the most promising are experimental conditions in 

mm.. 

Experimental conditions in other studies are rather 
close to the limiting case opposite to ([55]) . that is the 
case of high temperatures. For this reason we now con- 
sider, as in [l2|], the other limiting case 



T 
i/ n 



> 1, T < T c . 



(63) 



10" J and F 3 = 1.5 • lO" 3 respectfully. 



As before we assume here that e p ~ vqtiq. Then from 
([63]) follows that uj p (no) < T. 

In this case the following limiting expression for the dis- 
tribution function f® (u) p ) can be used (see, in this regard 
(HI]) and [3) 



T 

/p( w p)~ — j e p ~v Q n , w p (n )<T. (64) 

D p 

In order to simplify further calculations in (|40[) and 
(1411) it is convenient to change integration variable pto 
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variable z, that is introduced in accordance with the for- 
mula 



As a result of the change of integration variable we obtain 
the following expressions: 



E X (s) = 1 

X i i - 



Fa 



(s 2 - 1) Jo 

s 



2v (z) 



In 



dz g (s, z) 
\-v{z) 



s — v (z) 



(65) 



where now we have introduced the notations 

, . V2(3z + 1) \z(s 2 + l) + 2s 2 + ll 

9{s,z) = * 

V^(z + l)(z + 2'f 

r- z + 1 

v(z) = V2 



(66) 



y/z + 2' 

and the quantities F and a are still defined by (|44| and 

From the condition e\ (s ) = (see (2.19), flHU), (|6"T)l . 
{55])), written when s (a) s» ± (1 + 5s (a)), |<5s| -C 1, in 
the form of 



5s = Fa 



dz- 



2v (z) 



In 



(3z + l)(2z + 3) 
/2z (z + 1) (z + 2) 
l + v(z) + 5s 



l-v(z) + 5s 



(67) 



5s (a) = s (a) — 1, 



follows that the equation (|6TJ) has a solution only if 
Fa -C 1 (but a 3> 1 sec (f6"3")0 . The analytical solution of 
this equation does not exist. However, numerical analy- 
sis reveals that the result of integration with respect to 
z in the right-hand side of formula (l67l) almost does not 
depend on 5s if \5s\ -C 1. In this regard, the solution of 
(|6"7) in the case of <C 1 can be written as follows: 



5s = -Fad 2 , Fa < 1, 
where the constant di is given by 



(68) 



d-2 — I dz 
Jo 



yl+2 

2V2(z + l) 



In 



(3z + l)(2z + 3) 
sfTz (z + 1) (z + 2) 2 
VIT2 + V2(z + 1) 



VF+2- V2(z + 1) 



3.97. 

(69) 



temperature 



FIG. 5: Analytical dependence of dimensionless speed of 
sound s = uuo (q, a)/uoq on temperature (dimensionless quan- 
tity a) in case of "high" temperatures, plotted according to 
expressions JHH), for Ft = 3 ■ 1CT 4 , F 2 = 1CT 3 and 
F 3 = 1.5 ■ 10" 3 respectfully. 



Analytical dependence s (a) in case of o > 1 in ac- 
cordance with (|68p and (|69|) is shown on Fig. 4 for three 
different values of F. 

The value b(s) required in accordance with (I5T]) . (|52")) 
to calculate the damping rate of sound at high tempera- 
tures (see ([63]) ) is given by: 



Fa 
so (a) 7o 



9s 

rfzg(so(a),z) 
Fa 



gp («) + 1 

s (a) 



v 2 (z) 



x 1 



so (a) Jo 
so (a) 
2w (z, a) 



rfz 



s o ( a ) - v 2 ( z > a ) 
<9g (s (a) , z) 



(70) 



In 



<9s (a) 
s (a) + u (z, a) 



so (a) — v (z, a) 



where the functions g (s, z)and i; (z) are still defined by 
(|66| , and so (a) is given by (|67|) taking into account (|68|) . 
One can verify that the main contribution to b (s) the 
first term only yields. 

The damping rate 7 g of sound in BEC, calculated ac- 
cording to formulae (|51|) using (|70j) . in the main approx- 
imation with respect to Fa (see (|6"8"|) ) indicates a linear 
dependence on temperature: 



7 9 



w (?) 



u q 



7T (-7T + 6) 



Fa, 



(71) 



and that is a well known result (see for example [2, 12]). 
Deriving (jTTj) the following integral was used 



dz- 



(3z + 1) (2z + 3) 
/i(z + l) 2 (z + 2) 



7T + 6 



3/2 
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FIG. 6: Dependence of dimensionless damping rate 7 g /cJo 
on the temperature (dimensionless quantity o) in the case of 
'high' temperatures. The solid line shows the dependence (|71|) 
and the dashed one reveals the result of numerical computing 
based on the formulae (|50l) - (|52[) . Both lines were plotted for 

Formula ([71]) using (05]), (05) can be written as 

7? = ^~^-a sc Tq. (72) 

It should be mentioned that this expression differs 
prima facie from that given in 0, [HJ , where 

lq = ^-a sc Tq. (73) 

However, it is easy to verify that the numerical co- 
efficients in (IT2"j) and (173")) coincide with two significant 
digits. 

Fig. 6 shows the dependence of the dimensionless 
damping rate "fq/ujQ (q) at high temperatures region (a = 
(T/v n ) > 1 and T < T c ) plotted for F = 1CT 3 . It 
is apparent that the graph plotted according to formula 
(|TTj) or (|T2"j) (solid line) almost coincides with the graph, 
obtained by numerical computing based on formulae ([50")) 
- (152)) (dashed line). As it follows from the plot at Fig. 
4, the condition (-f q /ido) <C 1 of the existence of sound in 
dilute gas with BEC is well satisfied up the region where 
a 30. Naturally one should be confident that in a par- 
ticular physical system the condition T < T c is satisfied 
too in this range of a. 

V. CONCLUSION 

Thus we have reported the main results relating to the 
construction of microscopic theory of sound attenuation 
due to Landau mechanism in dilute gases with Bose con- 
densate. The analytical expressions of the propagation 
velocity and damping rate for sound in dilute gases with 
Bose condensate in the limiting cases of high and low 



temperatures were obtained. Shown that at high tem- 
peratures these expressions coincide with those ones ob- 
tained previously by other authors in various phenomeno- 
logical approaches. In the region of low temperatures 
the behavior of collisionless decrement of sound, obtained 
in the present paper, significantly differs from the same 
decrement, obtained by other authors. In our opinion the 
distinction is caused by our use of the evolution equations 
of investigated system that were obtained in the micro- 
scopic approach through the first principles. 

We now make a significant remark. At first glance, the 
presence of the general equations (14"5)) - (J5H) makes it pos- 
sible to find the parameters of propagation and attenua- 
tion of sound in present system in any temperature range 
(i.e. for all values a), at least numerically. Meanwhile, 
Fig. 3 and Fig. 4 do not display data for the 'intermedi- 
ate' temperature range, that is from a — 0.7 to a = 1.4. 
It is caused by the following circumstance. We have al- 
ready mentioned that in this range of values of athere are 
no analytic methods for solving dispersion equations (|4T))) 
- (I52|) . But as it turned out, in this temperature range the 
numerical methods, at least those we have used, also be- 
come uncontrollable. Hence one cannot trust its results. 
This is supposed to be due the fact that in this interval 
of a the point a « 0.847 is situated. Recall that it is 
the point where the value of sq (a) = ujq (q,a) /uoq (see 
(|49[) ) is equal to unity. As was previously mentioned the 
same value s gains also at the point a = (zero sound). 
It is that very neighborhood of a ~ 0.847 that gives the 
huge oscillations of sound damping rate in the outcome 
of numerical computation with insignificant changes of 
value a. However, one cannot be sure that these oscil- 
lations of damping rate reflect a real physical picture. 
The reason for doubt in outcome is in mentioned uncer- 
tainty in this temperature range of the numerical meth- 
ods used in the present paper. Numerical calculations 
are poorly controlled because of the slow convergence of 
the integrals in (|50|) - (f52j) due to the manifestation of 
non-analytic dependence on temperature of the disper- 
sion characteristics of the system in the neighborhood 
of point a w 0.847. Apparently more detailed study of 
the behavior of decrement (or maybe even increment!) 
of sound in the neighborhood of this point requires the 
use of more sophisticated numerical methods. Authors 
intend to address this issue soon. 

In this regard, the results in [32l are also of inter- 
est, where propagation and absorption of the transverse 
breathing mode of an elongated BEC were investigated 
experimentally in s7 Rb vapour. The authors of men- 
tioned study attempting to measure the damping rate 
of these modes at the temperature approximately equal 
40-60 nK had found that the behavior of the perturba- 
tion amplitude in this temperature region differs signifi- 
cantly from the behavior of the amplitude of a damped 
sinusoidal signal. The presence of such a phenomenon 
was suggested to be explained due to nonlinear effects 
in the propagation of studied modes in the system. We 
also would like to draw attention to the fact that the 
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FIG. 7: Dependence of dimensionless damping rate 7 9 /cjo on 
temperature (dimensionless quantity a). The solid line dis- 
plays the result of numerical computation for F = 3 ■ 10~ 4 
based on expressions ()50[) - (|52|) . Diamonds reproduce experi- 
mental points taken from [32] 



temperature range 40-60 nK corresponds to the values 
of quantity a in terms of the present article that are in 



the 0.7 — 1.0 range. The value of a is naturally calcu- 
lated according to formula (j44|) and using the values of 
the physical characteristics of the system [32|. In other 
words in the mentioned case [32| deals with the nearest 
neighborhood of the point a ss 0.847 where nonanalytic 
dependence of the dispersion characteristics on tempera- 
ture reveals. 

In conclusion we note a good correlation between our 
theory in the case of high temperature region with the 
experimental data [32]. As can be seen from Fig. 7, the 
results of numerical calculations based on formulae (|50[) 
- ([52]l for F = 3.2 ■ 10~ 4 (seefl44}) are in satisfactory 
concordance with experimental data 32]. Experimental 
data reveals the same good fit with dependence given by 
formula ([TTj) . 

The authors are naturally aware that a direct compar- 
ison of the results of the present paper with the exper- 
imental data [HI can hardly be considered as entirely 
correct. At least, the reason is that the authors of [32[ 
have to deal with trapped BEC. However, such a compar- 
ison yet again demonstrates the validity of statement [l2[ 
that the linear dependence of the sound damping rate on 
temperature should occur in the case of trapped BEC. 
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